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ABSTRACT 

Zaritsky and Lin (1997) have recently presented observations which can 
be interpreted as evidence for an intervening stellar population along the 
line-of-sight to the LMC, at a distance of about 34 kpc. The evidence is based 
on detected populations of red-clump and lower main- sequence stars which 
appear brighter than the principal distributions by about 0.9 magnitude, and 
which have similar fractional surface densities of ~6%. ZL97 suggest that the 
intervening population may belong to a low surface-brightness dwarf galaxy, or 
to the tidal debris tail of a disrupted galaxy. They derive an optical depth for 
gravitational lensing through the intervening system which is consistent with 
the observed value, thus suggesting that MACHOs comprise at most a small 
fraction of the Galactic dark halo. 

Such interpretation of the microlensing events toward the LMC has major 
implications for the nature of the dark matter, and thus should be tested in 
every conceivable way. We suggest a test which is based on the microlensing 
data themselves: if the observed events arise in the ZL97 intervening population, 
then all the lenses are at nearly the same distance, and have nearly the same 
transverse velocity. Thus, the predicted distribution of event durations should 
directly reflect the stellar mass function in the intervening population, which is 
probably not significantly different from that in the Galactic disk. 

We find that the first eight MACHO events toward the LMC are as consistent 
with a standard halo distribution of 0.5 M MACHOs (best fit), as they are 
with an intervening population having a transverse speed of ~ 240 km s _1 
and a Miller-Scalo stellar IMF. Monte Carlo simulations reveal that in order 
to completely rule out either the standard halo distribution of MACHOs or 
the ZL97 intervening population, about 40 events will be required if the IMF 
is rising below the H-burning mass limit, and ~ 140 if it is falling. If the 
intervening population hypothesis is correct then: (a) the currently observed 
range of event durations directly implies a lens mass range of ~ 0.06-0.8 M & , 
which brackets the mass of kO.2M where many investigations claim that the 
stellar mass function peaks, and (b) the event distribution provides valuable 
information on the low-mass end of the stellar IMF. 

Subject headings: gravitational lenses - dark matter - Galaxy: halo - Magellanic 
Clouds - stars: low-mass, brown dwarfs, mass function 
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1. INTRODUCTION 

Ground based observations of stars in the Large Magellanic Cloud (LMC) have recently 
revealed a population of red-clump stars that appear brighter by about 0.9 magnitude 
than the peak of the principal red clump distribution (Zaritsky & Lin 1997, hereafter 
ZL97). After considering and rejecting various possible explanations, ZL97 argue that these 
red-clump stars belong to a spatially localized intervening population of stars along the 
line-of-sight to the LMC, at a distance of about 34 kpc. ZL97 present corroborating evidence 
for this interpretation from HST observations of LMC stars (Holtzman et al. 1997). These 
data show a faint trace of a secondary lower main sequence population which is shifted 
by 0.9 mag and has a fractional surface density similar to that of the brighter red-clump 
stars (~6%). Using the derived distance and estimated projected density of the proposed 
intervening population, ZL97 derived an optical depth for gravitational lensing which is 
within a factor of two of the observed value, thus suggesting that MACHOs with masses in 
the range of ~ 0.1-1 M comprise at most a small fraction of the Galactic dark halo. 

These findings seem to support an earlier suggestion (Zhao 1997) that the observed 
gravitational microlensing events toward the LMC may be due to a low surface-brightness 
dwarf galaxy or the tidal debris of a disrupted system at the foreground of the LMC. Such 
interpretation of the microlensing data would have major implications for the nature of the 
dark matter, and thus should be tested in every conceivable way. Several potential tests 
have already been suggested: (a) looking for spatial correlations between the red-clump 
and the brighter red-clump stars in fields farther from the LMC center, as attempted 
by ZL97 who found suggestive support for the foreground population hypothesis; (b) 
analyzing variable stars in microlensing databases, as performed by Alcock et al. (1997b) 
who excluded the possibility of a dwarf galaxy along the l.o.s. to the LMC out to a distance 
of 30 kpc; examining surface-brightness maps within large areas around the LMC (Gould 
1997); and (d) looking for differences in the radial velocity distributions of the principal 
red-clump and brighter red-clump stars, as suggested by ZL97. 

In the present paper, we suggest a method utilizing the observed microlensing 
events themselves to discriminate between the possibilities that they arise in a standard 
halo distribution of MACHOs, or in an intervening stellar population at ~ 34 kpc. We 
demonstrate that the distribution of event durations should be different in these two cases, 
and examine the likelihood that the observed events are drawn from either population (§2). 
Using Monte Carlo simulations, we estimate the number of events required to safely rule 
out either hypothesis. In §3 we discuss possible implications of the results to the nature of 
the lensing population, and to the low-mass end of the stellar IMF. 
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2. DISTRIBUTION OF MICROLENSING EVENT DURATIONS 

The characteristic timescale of a gravitational lensing event towards the LMC is given 

by 

, _ 2r E _ 4[Gmx{l-x)d LMC ] 1/2 m 



where r E is the Einstein radius, x = dj d^uci d is the distance to a lensing object of mass m, 
dhMC is the LMC distance (~50 kpc; Feast & Walker 1987), and v t is the transverse velocity 
of the Einstein ring relative to the line-of-sight to the source. Therefore, the shape of the 
distribution of event durations depends on the density distribution of lenses along the l.o.s., 
p(x), their velocity distribution, and their mass function. These ingredient distributions are 
very different in the two cases between which we wish to discriminate. In the case of a dark 
halo distribution of MACHOs, p(x) is reasonably well determined, and the distribution of 
transverse velocities can be derived by assuming an isothermal velocity distribution, taking 
into account the Solar and LMC motions. The MACHO mass function is a priori unknown. 

By contrast, in the case of the intervening population proposed by ZL97, all the lenses 
are at nearly the same distance, di nt ~ 34 kpc, and have similar transverse velocities. The 
reason for the latter is the following: the angular size of the intervening system must not be 
much smaller than that of the LMC since otherwise microlensed stars would have appeared 
clustered in a limited area of the LMC (Maoz 1994). In addition, the mean surface mass 
density of the intervening population should not exceed 

47rG d int {d lmc - dint) V2.9xl0 V 

in order to be consistent with the inferred optical depth for microlensing toward the LMC 
(Alcock et al. 1997; hereafter A97). If the intervening system is bound, the dispersion in 
transverse velocities is cr int:± ~ (GM int j 'R ini ) x l 2 ~ (7rG f E int .R int ) 1 / 2 , where R int is the system's 
characteristic radius. This yields a inti ± ^24 (.R^/kpc) 1 / 2 km s^ 1 , which is very likely to 
to be negligible relative to the motion of the entire intervening system. In the case of an 
unbound or a tidally disrupted system, the total solid angle of the intervening structure must 
be much smaller than 4ir since otherwise its total mass would be ~4vrd 2 nt Si„i = 6xlO n M , 
which exceeds the entire mass of the dark halo. Thus, if the intervening population is to be 
relatively confined in space, its velocity dispersion must be much lower than the circular 
velocity in the halo, v c ~ 220 km s _1 (any possible large-scale velocity gradient along a 
tidal tail would not introduce much of a dispersion in the velocities of those intervening 
stars within the small solid angle of the LMC). We conclude that the distribution of event 
durations should directly reflect the mass function of the intervening population. Again, 
in contrast to the a priori completely unknown MACHO mass function, it is reasonable 
to expect that the stellar mass function in the intervening population is not significantly 
different from that in the Galactic disk. 
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2.1. Localized Intervening Stellar Population 

In the case of the intervening stellar population proposed by ZL97, we can approximate 
the lens distribution along the l.o.s to the LMC by 5(x — Xi nt ), where Xi nt — 0.68, and the 
transverse velocity distribution of the stars by 8(v — v t ), where 

v t = |v intj _L - v los (x int )\ , (3) 

v mt,_i_ is the transverse velocity vector of the intervening system, and vi os (x) is the transverse 
velocity of the line-of-sight to the source, at the distance of a lens, given by 

vios(x) = (l-x)v 0>1 + xv LMC ,± , (4) 

where v Qi j_ and vlmc,±. are the transverse velocity vectors of the Sun and LMC, respectively 
In a Cartesian coordinate system where v x , v y and v z are in the directions of the Galactic 
center, the Solar motion, and the Galactic North pole, respectively, the LMC motion with 
respect to the Galaxy is "Vlmc — (60 ±59, —155 ±25, 144 ±51) km s _1 (Jones, Klemola & Lin 
1994), the Solar motion with respect to the Galaxy is v Q = — v LSR +v c = (9, 231, 16) km s _1 , 
and the unit vector in the LMC direction is n LMC = (0.146, —0.826, —0.545). 
The differential event rate is 

dT = 2^(m)m l/2 v t R E (M Q ,x int )dm, (5) 

where rh = m/M & , Re is defined in Eq. (JJ), and \l/(m) is the present day mass function of 
the intervening stellar population, normalized such that / ^>(fh)mdfh = Y J i nt . Although the 
universality of the stellar mass function is yet to be proven, it is reasonable to assume that 
ty(m) is not significantly different from the stellar mass function in the Galactic disk. 

The understanding of the stellar mass function has advanced substantially since 
the concept of the IMF was introduced by Salpeter (1955), but at the low-mass end 
(~ 0.08-0.2 Mq), the results of various workers have been discrepant (e.g. see Miller & 
Scalo 1979; Simons & Becklin 1992; Kroupa, Tout & Gilmore 1993; Comeron et al. 1993; 
Tinney 1993). On one hand, there are claims for a turnover in the mass function at 
~0.2M Q , supported by observations of star forming regions (e.g. Hillenbrand 1997) and 
by possible evidence for a falling IMF at masses below the H-burning limit (~0.08M Q ) 
from observations of proto-brown dwarf clumps (Pound & Blitz 1995), thus suggesting that 
brown dwarfs exist in such small numbers as to make a negligible contribution to the mass 
density of the Galactic disk. At the other extreme, there are claims for a steeply rising mass 
function all the way down to the H-burning mass limit (e.g. Mera, Chabrier, & Baraffe 
1996), implying a substantial amount of brown dwarfs in the Galactic disk (see also Tinney 
1993). 

The shape of the low-mass end of the IMF is likely to be better determined in 
the next few years (e.g. due to deep IR searches for brown dwarfs), but at present we 



- 6- 



cannot rule out either possibility. Thus, we shall explore two mass functions which differ 
significantly at low masses: (a) the Miller and Scalo model (1979; hereafter MS79) where 
the IMF is fit by a lognormal distribution with a slope dlog^/dlogm = — (1 + logm), 
i.e. ty(m)dfh oc m urns over at O.1M ), and (b) the three-component power- 

law fit of Kroupa, Tout, and Gilmore (1993; hereafter KTG93), which is monotonically 
rising down to the H-burning limit, extended here down to a Jupiter mass 

3.5 rh- 13 0.001 < rh < 0.5 
#(m) oc { 1.9 rh- 2 - 2 0.5 < rh < 1.0 . (6) 

1.9 mr 2 - 1 1.0 < rh < oo 

The cutoff mass is introduced because the IMF cannot rise indefinitely, but the choice of a 
Jupiter mass is arbitrary. We note that the results will not be sensitive to the exact cutoff 
mass since the detection efficiency for very short events is low. 

Let us define t & as the crossing time of an Einstein ring due to one solar mass object 
at a distance of the intervening system, 

2r E (M Q ,x mt ) = / vt \ 1 

° v t V200 km s- 1 / y ' v ; 

where v t is defined in Eq. (|3|). Transforming dT /dm (Eq. ||) into dT/dt, we obtain for the 
MS79 mass function 

dT (1 \ -2pQg(£/t )+i] 

K (*AoJ > ( 8 ) 

and for the extended KTG93 mass function 

I 3.5 (t~A©)~°- 6 0.03 < i/t G < 0.707 
( '" oc { 1.9 (t/t )~ 2 - 4 0.707 <t/t e <l . (9) 



dt 



1.9 (tA©)" 3 ' 4 1 < */*© < 00 



Figure 1-a presents examples of both event distributions in the case of v t = 100 and 
300 km s _1 . 



2.2. Halo Distribution of MACHOs 

The event distribution in the case of a dark halo which consists of single mass MACHOs 
with a Maxwellian velocity distribution, taking the Solar and LMC motions into account, is 
(Alcock et al. 1996) 



dr l r 1 



oc 



dt mt 4 Jo 



dx p{x) r%{x, m) e -Q^™) e -M*V v <? I [P(x, m)] (10) 



where Q(x,m) = Ar%(x,m)/(f 2 vl), P(x,m) = Ar E ( X ,Tn)vi os (x) /(iv 2 ), I is the modified 
Bessel function, and vi os (x) = |vj os (a;)| is defined in Eq. M). 
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We shall consider here only a "standard halo" dark matter distribution of equal 
mass MACHOs (A97 present two examples of event distributions due to power-law mass 
functions, which closely resemble those due to equal mass MACHOs). Following Griest 
(1991), we take p(r)=p (a 2 +rQ)/(a 2 +r 2 ), where r is the Galactocentric radius, a = 5 kpc is 
the core radius, and ro = 8.5 kpc is the Sun's distance from the Galactic center. The density 
along the l.o.s. to the LMC is p(x) = p A' /(A' + Bx + x 2 ), where A'= (a 2 +r%)/d 2 LMC = 0.039, 
and B = -2(r /d LMC ) cos b LMC cos l LM c = -0.054. 



3. THE MACHO DATA vs. PREDICTED EVENT DISTRIBUTIONS 

Our ultimate goal is to rule out either the ZL97 intervening population hypothesis or 
the standard halo of MACHOs hypothesis by showing that it is extremely unlikely that the 
observed event durations were drawn from the predicted event distribution. In order to do 
that we need to take the detection efficiency of the relevant experiment into account. We 
shall test here the first eight MACHO events toward the LMC (A97) against the predicted 
event distributions, e(t) dT/dt, where e(t) is the detection efficiency for the first two years 
of MACHO data (Figure 2). The eight event durations are shown at the bottom of Figure 
1-b, along with examples of the predicted event distributions. 

Figure 3 shows the results of Kolmogorov-Smirnov tests (Press et al. 1992), where a 
is the probability of getting a K-S deviation worse than the observed value, namely, the 
significance level at which we can reject the null hypothesis that the data were drawn 
from a given distribution. For example, we can reject at the 0.1 significance level the 
hypotheses that the data were drawn from (a) the ZL97 intervening population with an 
extended KTG93 IMF when v t is either below 30 km s _1 or above 350 km s^ 1 , (b) the 
ZL97 intervening population with an MS79 IMF when v t <108 km s _1 , and (c) a standard 
halo distribution of single-mass MACHOs with m<O.15M . All remaining possibilities are 
currently consistent with the data, the more so for larger a. We see that the most likely 
MACHO mass in a standard halo model is 0.50 M Q , as concluded by A97, but the ZL97 
intervening population with an MS79 IMF and v t ^2A0 km s^ 1 is equally consistent with 
the present data. 

Eight events are clearly not enough for discriminating between the various hypotheses. 
The MACHO experiment has yielded so far an additional six events toward the LMC 
(Charles Alcock, private communication), but these are not yet publicly available. We 
performed Monte Carlo simulations to estimate the number of events required to safely 
rule out either the standard halo or the intervening population hypothesis, assuming the 
other hypothesis is correct (Figure 4). We find that a sample of ~40 events will enable 
to rule out a significant fraction of the parameter space with > 95% confidence level, but 
about 140 events will be required to completely rule out either hypothesis. However, this is 
the worst case scenario estimate for several reasons: (a) we assumed a detection efficiency 
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similar to that for the first two years of MACHO data, but e(t) is likely to improve with 
time, especially at long event durations due to the increasing baseline; (b) the simulations 
assume that the true event distributions are those with v t and m that best fit the present 
data, which need not necessarily be the case; (c) our understanding of the low-mass end of 
the stellar IMF may improve in the near future due to, for example, deep IR searches for 
brown dwarfs; and finally, (d) the permissible range of v t could be narrowed by subsequent 
investigations, as we discuss below. 

The transverse velocity of the ZL97 intervening population, Vi n t t j_, could generally 
be anything between zero and the escape velocity from the Galaxy, though not at equal 
probabilities: assuming a probability distribution for Vj n< similar to the velocity distribution 
function in an isothermal halo, P(vj nt ) d 3 Vi nt = (e - ^^) / (tc 3 ^ 2 v 3 ) d 3 Vi nt , we find that the 
probability distribution for v t (defined in Eq. is 

P(v t ) dv t = (v t /v c ) exp \-{v 2 t + (3 2 )/v 2 c ] I (v t /3/v 2 ) dv t , (11) 

where P = vi os (x int ) = 106 km s _1 . It peaks at v t = 236 km s _1 , but the distribution is quite 
broad. It may be possible to derive the transverse velocity of the intervening population by 
a careful re-examination of the data which led to the determination of the LMC's proper 
motion (Jones et al. 1994). Alternatively, if the radial velocity of the proposed intervening 
population is measured, then dynamical models may well constrain its transverse motion. 
Theoretical models may also yield predictions for the transverse velocity. For example, 
Zhao (1997) suggested that the lensing population toward the LMC may be a part of a 
tidal debris tail in the Magellanic great circle, in which case the tidal tail and the LMC 
are expected to have similar specific angular momentum. Neglecting the offset of the Sun 
from the Galactic center and substituting v int = x^t v lmc i n Eq. , we obtain for that 
scenario vt~201 km s" 1 . Therefore, the statistical test proposed in this paper may rule out 
various theories for the origin and nature of the intervening population with fewer events 
than estimated above. 



4. DISCUSSION 

If the observed microlensing events toward the LMC are due to the intervening stellar 
population proposed by ZL79, then the currently observed range of event durations (39-143 
days) directly implies a lens mass range of 0.06-0.8 (vt/200 km s -1 ) 2 M , which happens 
to bracket the mass of ~ 0.2 M & where many investigations claim that the stellar IMF 
peaks (see references in §2.1). On the other hand, if the observed events are due to a 
halo distribution of MACHOs, then the typical MACHO mass is ~O.5M (A97). These 
objects cannot be H-burning stars due to observational constraints from deep star counts 
(Gould, Bahcall & Flynn 1997), thus leaving the alternatives of primordial white dwarfs 
and sub-stellar mass black holes. It is intriguing that the ZL97 intervening population 
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hypothesis could provide a more "natural explanation" for the observed range of event 
durations. 

We propose a way to test not only the existence of a stellar population at the foreground 
of the LMC, but more importantly, the possibility that the observed microlensing events are 
primarily due to such an intervening population. Since the test relies on the shape of the 
stellar IMF, it could provide invaluable insight into the mass function below the H-burning 
mass limit if the ZL97 intervening population hypothesis is independently confirmed. Figure 
3 could already tell us that if v t ^,230 km s" 1 then a Miller-Scalo type of an IMF is much 
more likely than a monotonically rising one, and vice- versa if 1^130 km s -1 . 

Alcock etal (1996b, hereafter A96) placed strict limits on the contribution of 
low-mass MACHOs to the dark halo, based on the absence of detected events with 
durations of 0.1 ^20 days. If the observed events are due to the ZL97 intervening 
population, this null-detection provides an upper limit to the mass fraction in low-mass 
brown dwarfs. We can estimate the number of expected events due to, for example, 
10~ 4 ^ m <^ 10~ 2 (f t /200 km s -1 ) M objects by scaling the A96 results: denoting the 
fractional contribution of such objects to the total intervening mass by /, the expected 
number of short events is N exp ~ 35 / (t {, s / Thaio) , where 35 is the expected number of events 
due to MACHOs in this mass range, assuming they make the entire mass of a standard dark 
halo (A96, fig. 10), r obs is the microlensing optical depth inferred from the observations, and 
r haio is the predicted optical depth through a standard dark halo. Since T obs ~0.5r ha i o (A97), 
we expect N exp ~17f short events in the case of the intervening population hypothesis. For 
N o bs = 0, the Poisson 95% confidence level upper limit is 3 events, so we obtain an upper 
limit of / <z 18% to the mass contribution of 10~ 4 -10 -2 M objects to the stellar IMF. 
This limit will continuously improve with time as long as no short events are detected. 
Constraining the contribution of lower mass objects to the stellar IMF requires a detailed 
analysis which is beyond the scope of this paper since the detection efficiency of "spike" 
events does not depend only on i (A96). 

I thank Charles Alcock for discussion and comments, and Leo Blitz and Chris McKee 
for discussions. 
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Fig. 1. — (a) The predicted distribution of event durations toward the LMC for a standard 
halo model which consists of 0.5 M MACHOs (solid curve). The dotted (dashed) curves 
are for the intervening stellar population suggested by ZL97, assuming an MS79 (extended 
KTG) IMFs, respectively. Each set of curves corresponds to transverse velocities of v t = 100 
and 300 km s -1 . The triangles indicate the event duration due to a 0.08 M Q star (note that 
tocm 1 / 2 ). b) The predicted distributions for the same parameters, but after taking into 
account the detection efficiency for the first two years of MACHO data (Fig. 2). The points 
at the bottom show the durations of the first eight MACHO events (A97). All distributions 
are normalized to a unity area within 2 < t < 300 days, and take the Solar and LMC motions 
into account. 
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Fig. 2. — The detection efficiency for the first 2-year MACHO data (A97), kindly provided 
by Charles Alcock. 
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Fig. 3. — The results of K-S tests for the hypotheses that the first eight MACHO events 
were drawn from the predicted event distribution of (a) the ZL97 intervening population 
with transverse velocity v t and an MS79 IMF (dotted); same with an extended KTG93 IMF 
(dashed), and (b) a standard halo distribution of MACHOs of mass m. The lower a, the 
more inconsistent the data are with the predicted event distribution (see text). 



-14- 



■+-> 

CD 
> 
CD 



CD 



3 



-i-j 

CD 
> 
CD 



CD 



3 



100 



10 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

50 100 150 200 250 300 350 



V t [km/s] 



100 r 



10 ±* 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
MACHO mass [M ] 

Fig. 4. — Results of Monte Carlo simulations: (a) the mean number of events required 
to rule out (using a K-S test) a ZL97 intervening population with a transverse velocity v t , 
for an MS79 (dotted) and an extended KTG93 (dashed) IMFs, assuming the events arise 
in a standard halo distribution of 0.5 M MACHOs (the currently best fit to the MACHO 
mass). The upper and lower curves in each pair are for a = 0.05 and 0.1 significance levels, 
respectively, (b) the number of events required to rule out a standard halo of single-mass 
MACHOs of mass to, assuming the events arise in a ZL97 intervening population with (1) 
an MS79 IMF and best fit v t of 245 km s" 1 (dotted), and with (2) an extended KTG93 IMF 
and best fit v t of 140 km s _1 (dashed). See discussion in text. 



